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Abstract 


We study terminally attached self-avoiding walks and bridges on the simple cubic lattice, both by 
series analysis and Monte Carlo methods. We provide strong numerical evidence supporting a scaling 
relation between self-avoiding walks, bridges, and terminally attached self-avoiding walks, and posit that 
a corresponding amplitude ratio is a universal quantity. 


Keywords: self-avoiding walk; critical exponents; universal amplitude ratio; Monte Carlo; pivot algorithm 


1 Introduction 


A self-avoiding walk (SAW) on a lattice is an open, connected path on the lattice that does not revisit any 
previously visited vertex. Walks are considered distinct if they are not translates of one another. If we let the 
number of SAWs of n steps be Cn, it is known that lim„_>.oo n~^ logc,i = log p exists iflSTl . where p is the 
growth constant of self-avoiding walks on the lattice. This work will consider SAWs on the 3-dimensional 
simple cubic lattice Z^, with the vertices having integer coordinates for z = 0,1, • • • , n. 

The upper half-space H, is characterised by z > 0. An n-step bridge is a self-avoiding walk in the upper 
half-space that starts at the origin and is constrained so that ^(0) ^ ^{i) < ^(n) Q 

< i < n. We denote 

the number of n-step bridges starting at the origin by bn- It is known that lim^^oo log bn = log p, where 
p is unchanged from the corresponding value for SAWs lIT^ . 

A terminally attached self-avoiding walk, or TAW, is a SAW with one end anchored in the surface, but with 
the rest of the walk free in the upper half-space. Terminally attached self-avoiding walks are also referred 
to as half-space self-avoiding walks in the literature. Clearly TAWs are a superset of bridges, and a subset 
of SAWs, so have the same growth constant. The number of n-step TAWs starting at the origin is denoted 
by tn- 

The last subset of SAWs we wish to consider are arches, which are SAWs in the upper half-space with both 
the start- and end-point constrained to lie in the 2-dimensional surface z = 0. That is to say, = 0 = z^'^^ 
As the number of arches is bounded above by the number of SAWs and below by the number of self-avoiding 
polygons (SAPs), which are known to have the same growth constant as SAWs ll20l . it follows that arches 
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also have the same growth constant as SAWs. The number of n-step arches starting at the origin is denoted 
by an- 

The results on growth constants are the only results that have been proved. Nevertheless, it is universally 
accepted that the asymptotic behaviour of the above objects is given by: 


C-n 

~ A ■ 



(1) 

bn 

~ B ■ 



(2) 

in 

~ H 



(3) 


r^C ■ 



(4) 


for SAWs, bridges, TAWs and arches respectively. The existence of the various critical exponents has not 
been proved, but is universally accepted, and will be assumed hereinafter. 

All the above definitions also hold for the two-dimensional square lattice, where the upper half-space re¬ 
ferred to becomes the upper half-plane, and the originating surface is the line y = 0. 

Another critical exponent that needs to be defined is that which characterises the length of a SAW. All 
standard measures of the square of the length, such as the mean-squared end-to-end distance, the radius of 
gyration, the squared caliper span etc., behave as const, x where in two dimensions it is accepted that 
1 / = 3/4. This exponent can also be understood as the reciprocal of the fractal dimension, df. 

In two dimensions, the exponents are known exactly (assuming existence). In that case it is believed that 
7 = 43/32, 7i = 61/64, and 711 = —3/16. As far as we are aware, there are no published estimates for 7^, 
but as long ago as last millennium one of us (AJG) estimated the value of this exponent by series analysis to 
be 9/16. Subsequently, much longer series have been produced by Iwan Jensen, which enabled this estimate 
to be affirmed with much greater confidence. Somewhat later, Alberts and Madras (private communication) 
obtained the estimate 'jb = 9/16 for two-dimensional bridges through SLE arguments, subject to certain 
unproven assumptions, but this work was never published. 

In l|9l both SAWs spanning a strip and bridges were discussed, and comparisons made with conjectured re¬ 
sults from SLE 8 / 3 . The arguments given there allow one to predict 7;,, and they are reproduced, in summary, 
here. In |[28l it was explained why the measure of bridges starting at the origin and ending at x + iy should 
be y~^^^ f{x/y) for an (explicit) function / that decays exponentially. Summing over all values of x, it 
follows that the measure of bridges starting at the origin and ending at height h should be asymptotic to a 
constant times Then by integrating over height, the measure of bridges of height less than or equal to 

h should grow as 

Recall that for two-dimensional SAWs and bridges the exponent v = 1/d/ = 3/4. If we weight each bridge 
of length n by so that the normalised counts are given by it follows that normalised bridges of 

height h typically have steps. Hence the measure of bridges of n steps should be the same as the 

measure of bridges of height at most which gives for the measure. Equivalently, 

the number of bridges of exactly n steps grows like bn ~ const x yT' x = const, x //” x 

so that 76 = 9/16. 

These exponents are not all independent. There is a scaling relation, due to Barber |[T1, 

271 - 711 = 7 -h (5) 
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which holds independent of dimension. Very recently, Duplantier and Guttmann IfT^ have proposed the 
existence of another scaling relation giving the bridge exponent, 


lb = 711 + 


( 6 ) 


As expected, the exponents given above for two-dimensional SAWs satisfy both these relations. 

In three dimensions the most precise estimates we have are from the Monte Carlo work of one of us EHQ, 
with /i = 4.684039931 ± 0.000000027, 7 = 1.156957 ± 0.000009 and u = 0.587597 ± 0.000007, though 
the details of the Monte Carlo simulation for the estimation of 7 have not yet been published. For 71 there are 
a few estimates in the literature based on rather short series. Some 30 years ago, Guttmann and Torrie ifTTl 
estimated 71 = 0.676 ± 0.009, while the most recent Monte Carlo estimate is by Grassberger |[23l who 
estimated 71 = 0.6786 ± 0.0012. From Barber’s scaling relation this gives 711 = —0.3874 ± 0.0024. 


There is renewed interest in the values of the exponents 71 and jh in this case, as the former exponent arises 
in recent tests for the conformal invariance of 3d SAWs Il25ll26l . and the latter has been predicted by the 
recent scaling argument lITOll given above, and estimation of the exponent 7 ^ would be a useful test of that 
scaling relation. More precisely, Kennedy discusses the hitting density of a SAW in a connected domain 
when conformally mapped by a function / to another domain. The hitting densities with respect to arc length 
in the two domains are related by |/'(z)|^, where 6 is a critical exponent, identified by scaling arguments as 


b 


7 - 271 d 
2u ^ 2 ’ 

d 

2v^ 2' 


(V) 

( 8 ) 


Note that the above exponent estimates give the prediction h = 1.3296 ± 0.0020, compared to Kennedy’s 
direct Monte Carlo estimate b = 1.3303 ± 0.0003 lf2^ . For the bridge exponent the prediction from the 
pre-existing exponent estimates is ji, = 0.2002 ± 0.0024. 


In this paper we will estimate critical exponents and amplitudes associated with bridges and TAWs via series 
and Monte Carlo methods. 


1.1 Universal amplitude ratio 

Combining the two scaling relations, Q and Q, given in Sec.[T] we obtain 

271=7 + 76- (9) 

This exact relation suggests that the relationship between amplitudes corresponding to the LHS and RHS of 
the equation must be lattice independent, which in turn leads to a universal amplitude ratio. 

One interesting aspect of this ratio is that the objects involved are bridges and TAWs, which interact with a 
surface. Loosely speaking, the scaling relation Q may be interpreted as saying that the exponent associated 
with the number of pairs of TAWs of length n (which together have 2 free ends, and two ends anchored to 
a surface) is the same as the exponent associated with the number of pairs made up of a SAW and a bridge, 
each of length n. 
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In order for the lattice specific effects to cancel it is crucial that the boundary condition where the walk 
terminates on the surface is the same for the TAWs and both ends of the bridges. 

Now, the TAWs are anchored to the surface, and must satisfy the rule that the entire walk lies above this 
surface. In contrast, for convenience the bridges are defined fo be asymmefric, wifh an addifional sfep 
prepended fo one end fo liff if complefely ouf of fhe surface af one end. 

Thus we need modify fhe definifion of a bridge fo make if symmefric befween ifs fwo ends. The mosf nafural 
way fo do fhis is fo remove fhe frivial edge from fhe sfarf of fhe bridge. Convenienfly, a sfraighf-forward 
bijecfion exisfs befween modified bridges of lengfh n which satisfy ^(0) < < zi'^) for any 0 < i < n, 

and fhe usual bridges of length n + 1 which satisfy ^(0) ^ ^{l) ^ ^(n+l) fQj. Q 

< i < n + 1. Hence 

the number of modified bridges of length n is just bn+i, and this is the coefficient that must be used in the 
definition of the amplitude ratio. 


At the coefficient level, the above considerations lead to the following sequence of ratios 

( 10 ) 
( 11 ) 

where the amplitudes A, B, and H are from Q-([^, and AT is a universal constant. That is to say, we expect 
this ratio K to depend only on the spatial dimension of the lattice. 


Kn = 




Cn ■ bn+1 
K = lim Kr, = 


fxAB 

i /2 

fiAB 


n 


271 - 76-7 


(1 + 0 ( 1 ))) 


In terms of generating functions, the amplitude ratio is given as the limit 


lim 

X^Xc 


T{xf 

fi- C{x) ■ B(x) 


Universal constant. 


( 12 ) 


1.2 Confidence intervals 

We wish to highlight that the error estimates for the quantities calculated in this paper are not purely statis¬ 
tical, and hence there is a degree of subjectivity in obtaining them. In the analysis of relatively short series, 
which is the best that can be achieved for non-trivial 3d models such as SAWs, TAWs, and bridges, the 
interpretation of the results from series analysis can be quite subtle, and it is easy to be overly optimistic in 
interpreting the convergence (or lack thereof) of a sequence of improved estimates. For our Monte Carlo 
results we are approaching the large n limit, and consequently there is a lesser (but non-zero) amount of 
subjectivity involved. 

We have included plots where relevant to allow readers to judge for themselves how reliable our confidence 
infervals are. Loosely, one may think of the confidence infervals given in fhis paper as being roughly 
equivalenf fo a single sfandard deviafion. 
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1.3 Outline 


In Sec. we describe the generation of the series data for simple cubic lattice TAWs and bridges, which is 
then analysed in the Sec. Thereafter we describe the Monte Carlo computer experiment in Sec. and 
analyse the resulting data in Sec.[^ Finally, we discuss possible extensions and conclude in Sec.[^ 


2 Generation of series data 


We directly enumerated bridges and TAWs via a straight-forward optimized backtracking algorithm, whose 
running time was roughly proportional to the number of items being computed. 

To enumerate TAWs up to a certain length a recursive function was used that took the set of sites currently 
visited, the coordinates of the current site, and the maximum number of bonds left. The function first adds 
one to the tally of TAWs of the current length. It then, assuming there is at least one bond left, chooses each 
unused adjacent site (not allowing negative z coordinates) and calls the function again with that site as the 
current site, and one fewer sites remaining. 

This algorithm clearly has the function called exactly once for each TAW enumerated. As each iteration of 
the recursive function contains a small amount of code, with only one loop over adjacent sites (maximum 
6 ), its execution time is exactly proportional to the number of TAWs enumerated. 

Enumerating bridges uses a very similar algorithm. In this case the recursive algorithm keeps track of an 
extra variable; the maximum z value reached (^max)- Now instead of adding one to the TAW tally on each 
function call, a check is performed to see if the z coordinate of the current site equals Z ma x. 


As described, this algorithm would have execution time proportional to the number of TAWs, as exactly the 
same function calls would be made as for the TAWs. The performance can be improved to be proportional to 
the (significantly smaller) number of bridges by pruning the search tree when the current TAW cannot ever 
result in a valid bridge. This can be prohibitively expensive to compute exactly; a fast and very effective 
heuristic is to prune any branches when the 2 ; coordinate of the current point plus the number of remaining 
bonds is less than Zmax- This is somewhat conservative; it will not prune every dead branch as there might 
be an occupied site in the way; however in practice it works very well resulting in a runtime basically 
proportional to the number of bridges. 


We can reduce the work required to determine bn by enumerating irreducible bridges instead, which are 
the set of bridges that cannot be decomposed into a sequence of more than one bridge. Irreducible bridges 
can also be characterised as the set of bridges which have at least 3 bonds between any two layers, other 
than the first layer, which must have precisely one bond. Each bridge can in fact be uniquely decomposed 
into a non-empty sequence of irreducible bridges, and so the generating function for bridges, B{x), may be 
expressed in terms of the generating function for irreducible bridges, I{x), as 


B{x) 


i-i{xy 


( 13 ) 


Using this relation one may obtain bn from the series for the less numerous irreducible bridges. 
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Enumerating irreducible bridges requires a slight modification to the bridges algorithm. In this case yet 
another variable Zirr is kept track of in the recursive function: Zjrr is the largest z value below which the 
TAW is irreducible (it has at least three bonds in the 2: direction at each layer). To see if it is a valid irreducible 
bridge one now checks that the current 2; value, Zmax, and are all equal. is easy to keep track of; 
Zirr is incremented whenever a vertical step is made upwards from 2; = Z\„, as long as Z\„ < Z ma x. 

As described, this algorithm would have the same execution time as the bridges algorithm. Again, pruning 
can improve the speed to close to the significantly smaller number of irreducible bridges. The following 
conservative heuristic is used: 

• If 2 ^ ^irr ^ ^max? pfunc if the number of remaining steps is less than Z^ax — 2 : (this is the same as 
the bridges pruning); 

• If Zirr < 2 < Zmax? prune if the number of remaining steps is less than Zmax — '2 + l + 2(2 — Zirr) (this 
is due to the need to go back down to 2 = Zjrr, take one step sideways, and then return to 2 = Zmax)- 

We hnd that his heuristic works very well in practice, and so the running time is roughly proportional to the 
number of irreducible bridges. 

There were a variety of other optimizations used to give constant improvements, such as assigning each 
conceivably reachable lattice point an integer index and pre-computing adjacencies so that points were 
represented by integers rather than a tuple of integers. 

Symmetries were taken into account by enumerating all prefixes up fo lengfh roughly 8 (if varied somewhaf 
depending on whaf compufer fhe enumerafion was run on), canonicalizing fhem fo accounf for symmefries, 
and assigning fo equivalence classes. Then fhe recursive algorifhm described above was used on each 
equivalence class, and fhe resulfs mulfiplied by fhe number of members of fhe equivalence class. The 
enumerafions for each equivalence class were summed for fhe final resulf. 

If was fhen slraighl-forward fo make fhe algorifhm parallel (fo run on many differenl processors simulfa- 
neously). Each equivalence class can be compufed independenfly, and so were parcelled ouf fo mulfiple 
processors, wifh fhe only inferacfion being fhe final summation. 

The algorithm as described above is in no way specific fo fhree-dimensional sysfems, and indeed was made 
fo also work in ofher dimensions. Differenl lallice lypes could also be used, buf fhis may affecl fhe heurisfics 
used for pruning dead branches of fhe search free. 

The memory use is insignificanl; execufion lime is fhe conslrainl. 


3 Analysis of series 

As usual, we assume lhal fhe crilical behaviour of TAWs is given by 

r(x) = A(l-/i-x)-^i(l + A'(l-;U-x)^0, (14) 

n>0 
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and that of bridges by 


B{x) = ''^bn ■ ^ C{1 — fj, ■ x) '^'’{1 + C'{1 — ^ ■ x)^^). (15) 

n>0 

Here ( 6 „) is the number of n-step TAWs (bridges), counted up to translations, /r is the growth constant, 
for which we assume Clisby’s [i6i| estimate ^ = 4.684039931 ± 0.000000027. Here Ai = 0.528 ± 0.012 [4| 
is the leading correction-to-scaling exponent. A, A', C and C' are critical amplitudes (N.B., these are simply 
related to, but not the same as, the amplitudes for the series defined previously). 

In analysing the data, we have of course assumed the existence of the various critical exponents (which has 
not been proved), and have also used the best available estimate of the growth constant, given above, and the 
estimate of the correction-to-scaling exponent Ai. There are a number of estimates in the literature for this 
exponent, obtained by a variety of methods, including field fheorefical mefhods, series mefhods and Monfe 
Carlo mefhods. Esfimafes range from 0.47 up fo 0.57, wifh non-overlapping error bars. The recenf Monfe 
Carlo work of one of us |4| gives Ai = 0.528 ± 0.012, and if is fhis value fhaf we will use in our Monfe 
Carlo analysis. However our series work is less precise fhan fhe Monfe Carlo work, and for fhaf purpose 
faking Ai = 1/2 is sufficienfly precise. 

In Appendix]^ we give in Tablefhe series for TAWs up fo lengfh 26, and for bridges up fo lengfh 28. The 
bridge series is longer, as we acfually counfed fhe far less numerous irreducible bridges, and reconsfrucfed 
fhe bridge generating funcfion from fhe irreducible bridge generafing funcfion. We also recorded fhe heighfs 
of fhe irreducible bridges, and fhe full fwo-variable generafing funcfion for irreducible bridges, counting 
bofh lengfh and heighf, is given in Table [^in fhe appendix. The relationship befween bridges and irreducible 
bridges is given in ( fTS] ). 

We firsf analysed bofh fhe TAW and bridge series by using fhe sfandard mefhod of differenlial approxi- 
manfs lfT4ll . We inifially affempfed unbiased analyses of bofh fhe TAW series and fhe bridge series. The 
resulfs were disappoinfingly imprecise, despife choosing fo use fhird order differential approximanfs, which 
should, in principle, be able fo accommodafe fhe expecfed confluenl singularify strucfure. 

We had always infended fo use biased differenfial approximanfs, as fhese are expecfed fo give much more 
precise esfimafes of fhe exponenfs. Buf fhe disappoinfing lack of well-converged esfimafes in fhe unbiased 
case warns us nof fo be foo hopeful of obfaining excellenf resulfs in fhe biased case. Despife fhis, fhe 
esfimafes appeared tolerably well converged, and for TAWs we esfimafed 71 = 0.683 ± 0.005, while for 
bridges we estimated 7 ^ = 0.199 ± 0.004. The estimate of 71 jusf overlaps fhe earlier Monte Carlo estimate 
of Grassberger (231, and also jusf overlaps our much more precise Monfe Carlo esfimafe quofed below. 
We do nof know why fhis esfimafe of 71 is nof more accurafe. If is possible fhaf fhe amplifude of fhe 
correcfion-fo-scaling exponenf is quife large, buf we were unable fo confirm fhis from fhe available dafa. The 
corresponding esfimafe of fhe bridge exponenf % is in excellenf agreemenf wifh our Monfe Carlo esfimafe, 
given below, fhough fhaf esfimafe is again significanfly more precise. We do nof give more defails or fables 
of dafa here, as in fhe nexf subsection we give a more precise analysis. 

A completely differenl mefhod of analysis was also fried. We firsf calculafed fhe ratios of alfemafe co- 
efficienfs, where we have used alfemafe coefficienfs fo minimise fhe effecf of fhe anfi- 

ferromagnefic singularify locafed af x = —1/^. Then by sfandard ratio mefhod techniques l[T4l . one expecfs 
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the ratios to behave as 


~ /i IH-h 


,1+Ai + ^2 + 


(16) 


Here, for TAWs 5 = 71 — 1, and c and d are constants, while for bridges g = '^b — 1. Using the estimate of 
/i given above, we can then estimate the exponent g as the limit of the sequence gn defined by 


, \ const. const. 

9n= ( -l n~c/H-^ d-h 

g J n 


(17) 


We used the inbuilt fitting function of Maple to fit the elements of the sequence gn from n = rimin to 
n = Umax, where nmax = 28, to 


const. const. 

Qn = g -\ -Ai-1- 

n 


(18) 


with Umin Varying from 10 up to 21. The estimates of gn were reasonably well converged, and in this way 
we estimated g = —0.323 ± 0.003, or 71 = 0.677 ± 0.003. A similar analysis for the bridge series was less 
well converged. Fitting just to 


\ const. 

9n= { -l n~5rH-^ 

g ) n^i 


gave a decreasing sequence of estimates of Abridges, while fitting to 


const. const. 

gn = 9 -\ -Ad-1- 

n 


(19) 


( 20 ) 


as done for TAWs gave an increasing sequence of estimates of Abridges- Averaging corresponding entries 
gave a fairly stable sequence, from which we estimated Abridges = —0.802 ± 0.005, or 7 f, = 0.198 ± 0.005. 


Combining our results from the two methods of analysis, we give as our series estimates from this data, 
71 = 0.680 ± 0.003 and 7 f, = 0.1985 ± 0.004. Again, we do not give more details or tables of data here, as 
in the next subsection we give a more precise analysis. 


3.1 Series extension and subsequent analysis 

The analysis discussed above is based on the exact series coefficients given in Appendix [A| We can however 
obtain accurate approximations to the next five fo seven ferms of fhe series, which can fhen be used in our 
ratio analysis. Because every differenlial approximanf fhaf uses all fhe available series coefficienfs implicifly 
predicfs all subsequenf coefficienfs, we realised fhaf we could utilise fhis observafion fo calculafe, approx¬ 
imately, all subsequenf coefficienfs. Of course fhe number of significanf digifs decreases as fhe number of 
predicfed coefficienfs increases, buf as we show below, we can gel useful eslimales of fhe nexl six or seven 
coefficienfs. 

For every approximanf using all fhe known coefficienfs, we generaled fhe subsequenf six or seven coeffi- 
cienls. We observe fhaf fhe predicfed coefficienfs agree fo a cerlain number of significanf digifs among all 
the approximants, and we take these as our estimates. That is to say, assume we know the coefficients On 










for n G [0, A^max]- We then predict the coefficients aAfmax+ 2 ) • • • > «tVmax+7- Our estimate of each 

such coefficient is given by the average of the values predicted by the differential approximants. We reject 
obvious outliers, which are infrequent, and quote only those digits for which all approximants agree. So 
for example if the first 8 digits agree, and the coefficient is predicted to be a 20 digit integer, we will quote 
the coefficient as the 8 predicted digits followed by 12 zeros. Not surprisingly, we find fhe greafesf number 
of digifs are predicfed for wifh fhe number of digifs slowly decreasing as we generafed further 

coefficienfs. 

These predicfed coefficienfs are nof suifable for including in our differenfial approximanf analysis, as fhey 
were derived from lower order differenfial approximanfs. However for rafio fype analyses fhey are very well 
suifed, as discrepancies in say fhe sevenfh or eighfh significanf digifs will nof affecf fhe rafio analysis in fhe 
slighfesf. This is particularly useful in fhose sifuafions where we suspecf fhere mighf be a fuming poinf in 
the behaviour of ratios or their extrapolants with our exact coefficients, as these approximate coefficients are 
more than accurate enough to reveal such behaviour, if it is present. 

As a demonstration of this method, assume we only have 23 terms in the bridge generating function, and 
weTl predict the next 6 coefficients. In Table[T]we show the predicted and exact coefficients. (Note that the 
coefficient of is not exact, but is as predicted from a longer series, as discussed below). 


n 

Predicted bridges 

Actual bridges 

24 

7.115933 X 10^^ 

711593257794069 

25 

3.235637 x IQi'^ 

3235634079777801 

26 

1.472844 X 10^® 

14728414578753489 

27 

6.71105 X 10^6 

67110197685388181 

28 

3.06077 X 10^^ 

306074586987649389 

29 

1.39718 X 10^® 

1.397156394 x 10^® 


Table 1: Approximate, predicted coefficients of simple cubic lattice bridges of length n, compared to exact 
values. 


Note that in every case the predicted coefficient agrees with the exact coefficient to the order stated, with an 
uncertainty of a few parts in the last quoted digit. This demonstrates the utility of the method. 

In Table 1^ are the predicted coefficients from the full series, where, as explained above, we expect errors to 
be confined fo fhe lasf quofed non-zero digif. (Of course if is possible fhaf fhe lasf digif we are confidenf of 
is 0, buf nofhing is losf in glossing over fhis). 


We repeafed fhe basic rafio analysis described in fhe previous section, and show in Fig. [T] fhe esfimafors gn 
of fhe exponenf g for TAWs ploffed againsf Xj^/n^ as is appropriafe, see (17i. Sfrong oscillations can be 
seen, indicafing fhe presence of a significanf anfi-ferromagnefic ferm. 


We tiffed fhe dafa using fhe curve-tiffing feafures of Maple, simply by tiffing fhe dafa poinfs gn from 12 up fo 
Nma x — /c for fc = 0 ... 9, where Nmax = 35 in fhe case of bridges and is 34 in fhe case of TAWs, assuming 
gn is quadrafic in 1/ ^/n, as implied by (17 1 . The sequence elemenfs gn are shown in Table [^below. 


These sequences have been simply exfrapolafed, using every alfemafe ferm in fhe case of TAWs, fo give fhe 
limifs below. Nofe fhaf fhe exfra, approximafe, terms are (a) quite precise enough for this sort of analysis. 
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n 

TAWs 

Bridges 

27 

6.99275590 x 10^''' 

exaef 

28 

3.2418211 X 10^® 

exaef 

29 

1.5038283 X 10^^ 

1.397156394 x 10^® 

30 

6.9761657 x 10^^ 

6.38289508 x lO^® 

31 

3.237937 x lO^^ 

2.91825221 x 10^^ 

32 

1.502907 X 10^1 

1.33518243 x lO^o 

33 

6.979108 X 10^1 

6.1129712 X 10^0 

34 


2.8005385 x lO^^ 

35 


1.2837858 x 10^2 


Table 2: Approximate, predieted eoeffieients of simple eubie lattiee TAWs and bridges of length n. 




Figure 1: Exponent sequenee gn for TAWs 


and (b) move the data into a regime where the extrapolations are mueh elearer than those using only the 
exaet eoeffieients. The extrapolation method was based on the observation that the numerieal gaps between 
sueeessive (or alternate, as appropriate) eoeffieients apparently deerease as a geometrie sequenee, so this 
was Slimmed. 

As a result of this analysis we find 71 = 0.676 ± 0.002, and 7 ^ = 0.199 ± 0.002. These are eonsisfenf wifh, 
buf more preeise fhan, fhe resulfs obfained in fhe previous seefion, using only fhe exaefly known eoeffieienls. 
In subsequenf seefions we diseuss fhe Monfe Carlo resulfs, whieh are subsfanfially more preeise fhan fhe 
series resulfs. However, as we show below, fhere is one properly whieh is aeeessible from series analysis 
whieh is nol easy lo delermine by Monfe Carlo analysis, and fhaf is fhe behaviour of spanning bridges in a 
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Figure 2: Exponent sequence gn for bridges 


k 

9n TAWs 

gn bridges 

9 

-0.34067 

-0.86711 

8 

-0.34120 

-0.85167 

7 

-0.33687 

-0.83958 

6 

-0.33707 

-0.83011 

5 

-0.33429 

-0.82278 

4 

-0.33432 

-0.81715 

3 

-0.33244 

-0.81287 

2 

-0.33240 

-0.80970 

1 

-0.33207 

-0.80739 

0 

-0.33098 

-0.80577 

Fimit 

-0.324 ± 0.002 

-0.801 ±0.002 


Table 3: Sequences of exponent values gn for TAWs and bridges. 


slab of given thickness, at the critical temperature, as the slab thickness increases. We discuss this in Sec. 3.3 


below. 


3.2 Universal amplitude ratio 


From our series, and known series for 3d SAWs MM, we have calculated Kn for re < 27 using exact 
coefficients, and for re < 34 using the last seven approximate coefficients. From the asymptotic form of the 
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( 21 ) 


coefficients, we expect the asymptotic form of the amplitude to be 


Kn-^K 1 + 


h 


n‘- 


k2 

+ — + 


n 


Accordingly, we extrapolated Kn against 1/n'^i, using our approximate value Ai = 1/2. As can be seen 
in Fig. 1^ the sequence is monotone decreasing, as far as we have data, with the last entry being 
i ^34 = 0.61601..., which suggests that this is an upper bound on K. The ratio plots against Ijy/n ex¬ 
hibit significant curvature, consequently we could only form the estimate K > 0.609. Combining these 
numerical confidence bounds, our first estimate is K = 0.6125 ± 0.0035. We also extrapolated Kn against 
1/n, and those plots exhibited similar, though reduced curvature. The increased linearity when plotting 
against 1/n suggests that the amplitude of the leading correction-to-scaling term is rather weak, and that the 
sub-dominant behaviour for n < 34 is largely given by the 0(l/n) term. Accordingly, we eliminate this 
term by forming the estimators 

( 22 ) 

Note that alternate terms are used to reduce oscillations caused by the anti-ferromagnetic singularity present 
in all three generating functions at x = — l//r- We expect 


^'-(1 + ^17 +»(/))■ <«) 

We show in Fig.j^a plot of Kn against Ijy/n, and the importance of the extra estimated terms is now seen, 
as we are just getting the first clear indication of a turning-point in our estimates. If we plot the odd and 
even terms separately, this is even clearer. From this figure, and the separate odd-term and even-term plots 
(not shown), we can give the refined estimate K = 0.6140 ± 0.0020, which is shown on the graph. 



Figure 3: Amplitude estimators Kn and Kn for universal amplitude K. Our best estimate from the series 
analysis K = 0.6140 ± 0.0020 is also shown. 

A similar analysis for the square lattice was also conducted. There we have much longer series available, and 
the dominant correction is 0(l/n), as in two dimensions Ai = 1.5. In that case we obtained K 2 d ~ 0.57280 
where we expect the error to be confined to the last digit. 


12 






3.3 Behaviour as strip width grows. 


The data we have generated also allows us to estimate the exponent characterising the behaviour of the 
generating function of bridges, at the critical point, which span a strip of width T as T —)■ oo. For two- 
dimensional SAWs, under the assumption that these are describable in the scaling limit by SLE 8 / 3 , one 
has |[28]l 


Bt{Xc) ~ 


const. 

rV4 


For the simple cubic lattice, we expect similar behaviour, so that 


Bt{Xc) ~ 


const. 

j'd 


(24) 


(25) 


From the data for simple cubic bridges, given in the appendix, we can construct the first 28 terms in the 
generating function for bridges in a strip of width T. Note that such generating functions have as their first 
non-zero term, as one needs at least T -steps to span a strip of width T. Note too that the generating functions 
will diverge at a critical value xt > Xc, where Xc = l//r is the bulk critical value. So evaluating Bt{xc) 
means we are evaluating these ordinary generating functions at a value inside their radius of convergence. 


In order to do this as effectively as possible, we have used Fade approximants. For each width T, we have 
constructed a number of Fade approximants utilising all, or nearly all, of the coefficients, and evaluated 
these at Xc. As T increases, the spread among the approximants for a given value of T increases, and we 
found we could only get useful estimates for T <9. These were 0.6967886, 0.5565774, 0.4720324, 0.4144, 
0.3728, 0.3399, 0.3132, 0.2915, 0.273 for T = 1,..., 9. We do not quote errors, but expect the error in each 
case to be confined to the last quoted digit. A log-log plot of these data, which should have gradient —9, 
shows slight curvature. Accordingly, we calculate the local gradient given by 


Q _ log Bt{xc) - \ogBT-l{Xc) 
^ ~ log(r) - iog(r - 1 ) 


(26) 


We plot 6t against reflecting the correction-to-scaling exponent around 0.5, which gives a straight 

line (with minor oscillations). Extrapolating this with a straight-edge, which is an appropriate level of 
sophistication given the paucity of data we have, gives a best estimate of 0 = 0.7. If the data is representative 
of the asymptotic behaviour, as it appears to be, quoting a confidence interval of ± 0.1 is appropriately 
conservative. 


We can also give an heuristic argument, based on the ideas given in Beaton et al. Q, as to the expected 
value of this exponent in terms of known exponents. We first argue that the measure of bridges of length 
< n should scale as where 7 ;, is the bridge exponent introduced above. Now normalised bridges of 
height T typically have T^/'^ steps. So the measure of bridges of height < T behaves asymptotically as 
(j' 1 /!^) 76 . Then bridges of height exactly T should scale as ^ T“®, i.e. 

6 ' = !-—. (27) 

V 


In two dimensions v = 3/4, and 7 ;, = 9/16, so 0 = 1/4, as expected. In three dimensions we have 
V = 0.587597 ± 0.000007, and our series estimate for 7 ;, is 0.199 ± 0.002, so the above equation gives 
9 = 0.661 ± 0.004, which is consistent with the direct estimate. 
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4 Generation of Monte Carlo data 


For our computer experiments we utilise the pivot algorithm, a Markov chain Monte Carlo scheme which 
allows for the efficient sampling of self-avoiding walks and related objects. The pivot algorithm was in¬ 
vented by Lai ll27l but first studied in detail by Madras and Sokal |[^ . In particular, we use the recent 
implementation of one of us Hill, which improved on earlier important work by Kennedy Il2^ . 

The pivot algorithm is a method to sample self-avoiding walks and TAWs of fixed lengfh, which nafurally 
gives esfimafes for quanfifies associafed wifh fhe mean size and shape of a SAW, such as v. However, if is 
less obvious how fo exfend fhe mefhod fo calculate quanfifies which inherenfly depend on differenf lengfhs. 
We gef around this difficulty by calculating ratios between different kinds of walks. The basic method is 
to efficiently sample walks in some larger set, and then at each time step test whether the walk belongs to 
some smaller set. For the sake of efficiency if is crucial fhaf fhe rafio of fhe sizes of fhese fwo sefs is nol 
exponentially small; in pracfice fhe rafio is a negative power of fhe lengfh n. 


The fwo quanfifies we esfimafe are fhe probabilify fhaf a self-avoiding walk is a TAW, and fhe probabilify 
fhaf a TAW is a modified bridge which has had ifs firsf sfep removed. We choose fo esfimafe fhis rafio for 
fwo reasons: if is more sfraighf-forward for our implemenfafion of fhe pivof algorif hm fo adapf if in fhis 
wa>[^ and also because if facililales fhe esfimafion of fhe amplifude ratio K from Sec. 3.2 fhaf we anticipate 
is universal. 


For fhe firsf computer experimenf, we sampled self-avoiding walks via the pivot algorithm, and for each time 
step we tested whether the SAW is in fact also a TAW. Our observable is the indicator function Xh{^)^ which 
is 1 if cj is a TAW, and zero otherwise. When we take the expectation of Xh with respect to the uniform 
distribution on the set of SAWs of length re, from ([T]l and Q we have 

{Xk)n = 

^^nre7i-l 


Thus estimating the mean value of this observable for a variety of lengths allows us to estimate the difference 
between the exponent for TAWs, 71 , and the exponent for SAWs, 7 . 

Similarly, for our other computer experiment we sampled TAWs via the pivot algorithm (N.B., the pivot 
algorithm is ergodic in this case, even though no proof has appeared in the literature as far as we know), and 
for each time step we tested whether the TAW is indeed a bridge with its first step removed. Our observable 
is the indicator function Xb{^)^ which is 1 if w is a bridge with its first step removed, and zero otherwise. 
When we take the expectation of Xb with respect to the uniform distribution on the set of TAWs of length re 

'Technical aside: this is because we only need to check that the ends of the walk lie on opposite faces of the minimum bounding 
box, without needing to check if additional sites might lie on the faces. 


(28) 

(29) 

(30) 
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{'Hn), from Q and Q we have 


/ \ bn+l 

{XbJHn = - 


Fin 




(31) 

(32) 

(33) 


This allows us to estimate the difference between the bridge exponent jj, and the TAW exponent 71 . 


One key detail of these simulations, which is similar to the calculation of the connective constant p in Il 6 l. 
is that if the pivot site is selected uniformly at random then this would result in an integrated autocorrelation 
time for the ratio which is 0(n). The reason that this is true for Xh^ in brief, is that a positive fraction of 
SAWs must have at their end a configuration which is of some fixed lengfh k fhaf cannof ever be affached 
fo fhe surface. The shorfesf possible example for fhe beginning of a self-avoiding walk which can never 
be a TAW on fhe square laffice is shown in Fig. Whenever fhe beginning of a walk assumes such a 
configurafion, if musf necessarily be fhe case fhaf Xh{^) = 0 - 



Figure 4: Minimal walk beginning which can never be a TAW on fhe square laffice (solid line) wifh a 
possible exfension (dashed line). 


If we were fo sample from fhe n possible pivof sites uniformly af random, fhen fhe probabilify of selecting 
one of fhe firsl k sifes is k/n. Thus we expecf fhaf if musf fake on average 0(n) time sfeps before a pivof site 
is selecfed sufficienfly closely fo fhe end fo have any chance in changing Xh from 0 fo 1. A similar argumenf 
applies fo fhe sampling of TAWs and fhe esfimafion of Xb- 

For each of fhe computer experimenfs, insfead of sampling pivof sifes uniformly af random on fhe n-sfep 
SAWs and TAWs we preferentially sample close fo each of fhe ends, using fhe following fhree disfribufions 
for pivof sifes f G [0, n — 1], each selecfed wifh probabilify 1/3: 


• Infeger i sampled uniformly from fhe interval [0, n — 1]; 

• Real X sampled uniformly from [0, log(n + 1)), i = [e^J — 1; 

• Real X sampled uniformly from [0, log(n + 1)]), / = n — [e^J. 


This choice of pivof sife sampling disfribufion guaranfees fhaf all lengfh scales wifh respecf fo fhe disfance 
from each end are rapidly sampled. We sample uniformly from all symmefries of fhe simple cubic laffice, 
excluding fhe identify. 
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SAWs were initialised using the pseudo-dimerisation algorithm described in [5 ], followed by approximately 
20n successful pivots, with pivot sites chosen uniformly at random, to equilibriate the Markov chain. 

The TAWs were also initialised using pseudo-dimerisation, but an additional sequence of pivot moves was 
performed until the SAW was in the half-space. For further equilibration, we wanted to ensure that there 
was no bias due to the surface, and so instead of sampling pivot sites uniformly we sampled pivot sites 
i G [0, n — 1] with probability 1/2 from the following two distributions: 


• Integer i sampled uniformly from the interval [0, n — 1]; 

• Real X sampled uniformly from [0, log(n -|- 1)), i = [e^J — 1. 


We performed approximately 40n successful pivots in this case. 

Each computer experiment was run for 21 different values of re, with re = 2^ — 1 for /c = 9,10, • • • , 25, with 
additional data for re at re = 723,1447, 2895, 5791, as we found that there was more information available 
for our estimates from smaller values of re. Approximately 20000 CPU hours were used in each experiment, 
divided equally between the different lengths, on SunFire X4600M2 machines with 2.3GHz AMD Opteron 
CPUs. These data are reported in Appendix [B| 

In addition, we invested a relatively small amount of computer time gathering data for shorter lengths. 
We do not report these result here, but show them in graphical form below in Fig. to illustrate strong 
corrections-to-scaling. 


5 Analysis of Monte Carlo data 


For our analysis we used the method of direct fitting. We fitted our data from Appendix]^ with the asymp¬ 
totic forms for each of the series, using in each case the term with the leading correction-to-scaling involving 
Ai. We have also performed fits of the dominant term only, as an alternative, robust, but less accurate method 
for determining the values of critical parameters. These more robust fits confirm fhaf if is reasonable fo fif 
the Ai term in each case. 

We used weighted least squares to fit all values for re > remin^ and plotted the fits against an appropriate 
power of remin which was of the same order as the first neglected correction-to-scaling term. We only 
required linear fits, and used the statistical programming language R as our tool. We performed fits with 
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nmin < 16383, as statistical noise for larger rimin rendered these fits useless. The fitting forms we used are: 


{Xh)n — 
log(X/i)n ~ 
{xb)nn = 


^og{xb)n„ ~ 


{Xh)n 


{Xb)Hr, 
1 {Xh)n 


r\j 


t 


H 


— ~ ( 1 + 


const. 


n 


+ 0 - 


n 


/ M / e X 1 ^ const. 

(71 - 7 ) log n + 6h) + log — + , . 

A (n + 4)^i 


+ 0 


^ ~ ( 1 + 
tn. H 


const. 

n'^i 


+ 0( - 
n 


/ X, X r X , const. 

(7» - 7l) log{n + h) + log ^ 


+ 0 - 


n 




^^ 271 - 76-7 ( 1 + 


const. 

n^i 


+ 0( - 
n 


bn+lC-n 

(271 - 76 - 7 ) log n + log iT + + O ( - ) 

\n J 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


The neglected next-to-leading correction-to-scaling terms have powers —1, —2Ai, and —A 2 , which are all 
approximately —1 (which is why we write 0(l/n) as a shorthand). We chose not to attempt to fit these 
next-to-leading corrections, due to the competition between these terms which makes it difficult to sensibly 
interpret any fitting procedure. In addition, we neglected the anti-ferromagnetic singularity which we expect 
to be small for the values of n considered here. Fortunately we have the luxury of being able to study the 
large n regime where we expect all of these correction terms to be small compared to the leading 
correction. 


As was the case for the series analysis, we biased the leading correction-to-scaling exponent Ai. However, 
since our Monte Carlo estimates are more accurate, we used the best available estimate Ai = 0.528 ± 
0.012 14] rather than the approximate value Ai = 0.5. We performed fits for Ai = 0.516, 0.528,0.540, but 
only plotted the confidence infervals for fhe cenfral esfimafe. We did fhis fo gain an understanding of the 
relative contributions of statistical error and the systematic error due to biasing the fits with an imprecise 
value for Ai. 


The parameters 5h and 5h were introduced by hand to reduce the slope of the fits in Figs. [5]|^ We found 
that 6h = 2 and <^6 = 3 do this job adequately, which made it easier to extrapolate our estimates to n = 00 . 
Effectively, this perturbation provides a crude method of approximately cancelling out the next-to-leading 
corrections to scaling. This trick is not useful for the fits involving the universal amplitude ratio, as in this 
case the leading exponent is zero which means the corrections introduced have small amplitude. 


For the critical exponents corresponding to the amplitude ratio, we performed fits using (391 to verify that 
the combination of critical exponents is very nearly zero. For the amplitude ratio, we assumed that the 
combination of exponents in (1^ is exactly zero, and used ([38]) to directly estimate the amplitude. 


In each case we have confirmed fhaf fhe goodness-of-fif is approximafely 1 in fhe large n regime, which 
indicafes fhaf fhe fiffing form is sensible and fhaf sub-leading ferms are losf in fhe sfafisfical noise. 


In Figs and [6 we plof estimates of exponenfs and amplifudes from (351, obfaining esfimafes for 71—7 
and ^. In Figs ^ and we plof esfimafes of exponenfs and amplifudes from (37 1 , obfaining esfimafes for 


76 — 71 and ^. Note fhaf fhe scaling relation in 


9b implies fhaf 71 — 7 = 76 — 71 ■ 
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Figure 5: Estimates of 71 — 7 from fits to Monte 
Carlo estimates of ratios tnicn, with correction- 
to-scaling exponent Ai = 0.516,0.528, and 
0.54. Our extrapolated estimate, 71 — 7 = 
—0.479290(15), is shown in the plot. 




7 + 4 )-' 

Figure 6 : Estimates of log ^ from fits to Monte 
Carlo estimates of ratios tnIcn, with correction- 
to-scaling exponent Ai = 0.516,0.528, and 
0.54. Our extrapolated estimate, log ^ = 
0.32520(30), is shown in the plot. Hence ^ = 
1.38431(41). 
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Figure 7: Estimates of 7 ^ — 71 from fits 
to Monte Carlo estimates of ratios bn+iltn, 
with correction-to-scaling exponent Ai = 
0.516,0.528, and 0.54. Our extrapolated esti¬ 
mate, 7 fe — 71 = —0.479315(20), is shown in 
the plot. 
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Figure 8 : Estimates of ^B/H from fits 
to Monte Carlo estimates of ratios bn+i/tn, 
with correction-to-scaling exponent Ai = 
0.516,0.528, and 0.54. Our extrapolated esti¬ 
mate, log ^ = 0.81220(30), is shown in the 
plot. Hence ^ = 2.25286(68). 


A A, = 0.540 
o Ai =0.528 
V Ai =0.516 



In Fig.we plot the universal amplitude ratio as a function of n, to illustrate the grave difficulty in extracting 
reliable estimates from short series. In this case no fits are performed, we simply plot the amplitudes against 
the leading correction to scaling which corresponds to Ai. One can see that the sub-leading corrections to 
scaling are extremely large, and overwhelm the leading Ai corrections to scaling until n is of the order of 
500 or so. To be able to extract reliable estimates from the low order coefficients one must therefore take 
into account the next-to-leading corrections. 


In Fig. 10 we show our estimates of 271 — 7 b — 7 from fitting (391. The fits strongly suggest that this 
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Figure 9: Plot of Kn = (bn+iCn) for n > 20 to show the strong next-to-leading corrections to scaling. 
It is impossible to reliably extrapolate these data for n ^ 500 without taking the next-to-leading corrections 
into account. 


combination of exponents is indeed identically zero. In Fig. 11 we show our fits of (38 1 with 271 — 76 — 7 
biased to be zero, i.e. we assume that the scaling relation (0 is correct. By doing this we can obtain far more 
accurate fits of this putative universal amplitude ratio than we were able to for the other amplitude ratios. 
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V Ai = 0.516 
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Figure 10: Estimates of 271 — 76 — 7 from 
fits to Monte Carlo estimates of ratios Kn = 
with correction-to-sealing expo¬ 
nent Ai = 0.516,0.528, and 0.54. Our extrapo¬ 
lated estimate is 271 — 7 & — 7 = 0.000025(30). 


Figure 11: Plot of amplitude estimates from 
fits to Monte Carlo estimates of Kn = 
t^/(bn+iCn), biased to have 271 - 76 - 7 = 0 , 
and with correction-to-scaling exponent Ai = 
0.516,0.528, and 0.54. Our extrapolated esti¬ 
mate isK = 0.614730(30). 
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5.1 Summary of results 


We summarise our Monte Carlo results in Tableproviding comparisons between different estimates of 
271 — 7 b — 7 and K from fitting data from the two computer experiments separately, and together. Reas¬ 
suringly, we find in each case that the confidence intervals overlap. Note that the direct estimate of K is an 
order of magnitude more accurate than the combined estimate, due to the biasing of the exponents to satisfy 


Quantity Estimate 


71-7 

-0.479290(15) 

7 b - 7 i 

-0.479315(20) 

271 - 7 b - 7 

0.000025(25) 

271 - 7 b - 7 

0.000025(30) 

H 

A 

1.38431(41) 

fiB 

~Tr 

2.25286(68) 

K 

0.61447(26) 

K 

0.614730(30) 


Source 


Fig. 

Fig. 


Fig. 1^ and Fig. 
Fig. 


10 


Fig. 

Fig. 


Fig. and Fig. 
Fig. 11 


0 


0 


Table 4: Summary of estimates made via analysis of Monte Carlo data. 


In addition, we combine our direct Monte Carlo estimates with previous exponent estimates 7 = 1.156957± 
0.000009 ||71 and u = 0.587597 ± 0.000007 [4J. We then obtain estimates for the TAW exponent 71 = 
(71 “7)+ 7 = —0.479290(15) + 1.156957(9) = 0.677667(17), and the bridge exponent 7 b = (7b —7i) + 
(71 - 7 ) + 7 = -0.479315(20) - 0.479290(15) + 1.156957(9) = 0.198352(27). N.B., we combine the 
confidence intervals as if they were independent statistical estimates of a single standard deviation. 


The Monte Carlo estimates are consistent with the series analysis estimates of Sec. | 3 . 1 | ( 7 i = 0.676 ± 0.002 
and 7 b = 0.199 ± 0.002) but are dramatically more precise. 

From the scaling relation Q, we calculate 711 = 7 b —= 0.198352(27)— 0.587597(7) = —0.389245(28). 

From the scaling relation ([^, we calculate b = — 7 b/( 2 i/) + dj2 = 1.331218(23), which may be compared 
with the direct estimate due to Kennedy of 6 = 1.3303(3) Il26ll . 


For the exponent 6, the hand-waving argument in Sec. |3.3| gave us the scaling relation ( [27] ), from which 
we obtain 9 = 1 — 7b/j^ = 0.662435(46). This indirect estimate of 9 should be compared with the direct 


estimate of 0 = 0.7 ±0.1 and the indirect estimate of 0 = 0.661 ± 0.004 from the series analysis in Sec. 3.3 


We summarise these results in Table In each case, our estimates improve significantly on those previously 
available in the literature. 
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Quantity 

Estimate 

71 

0.677667(17) 

lb 

0.198352(27) 

111 

-0.389245(28) 

b 

1.331218(23) 

9 

0.662435(46) 


Table 5: Summary of estimates for critical exponents from Monte Carlo. 


6 Discussion and conclusion 

6.1 Possible improvements and extensions 

For the enumerations, we utilised a straight-forward backtracking algorithm despite the fact that some im¬ 
proved algorithms have been developed in recent years for self-avoiding walks for d = 3. The three innova¬ 
tions are the use of the lace expansion and the two-step method [8], and the length-doubling algorithm lISTl . 

Of these methods, it is unlikely that the lace expansion would help as this works by enumerating configura¬ 
tions like self-avoiding polygons instead of SAWs, and there would not be as much of a gain for TAWs and 
bridges. 

The length-doubling algorithm is extremely promising, but would be challenging to adapt to the present 
problem. If this could be done, then it may be possible to extend the series past the existing n = 36 for 
SAWs imi . 

Finally, the two-step method could be adapted quite easily and used to increase the length of the series for 
the same computational effort, perhaps by a few terms. 

Unfortunately it is unlikely that the extent of the improvement would be sufficient to be competitive with 
Monte Carlo, even in the case of the length-doubling algorithm. 

To improve the Monte Carlo results it would be possible to perform the same computer experiment with 
weakly self-avoiding walks (also known as the Domb-Joyce model) where the weight factor is chosen so 
that the leading correction-to-scaling term corresponding to Ai has small amplitude. 

In future we will perform further Monte Carlo computer experiments for a variety of two- and three- 
dimensional lattices to test our hypothesis that K is indeed a universal amplitude ratio. We have some 
preliminary numerical evidence that this ratio is the same for the simple cubic, face-centred cubic, and 
body-centred cubic lattices. 
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6.2 Conclusion 


We have obtained high precision estimates of a variety of exponents and amplitudes associated with termi¬ 
nally attached self-avoiding walks and bridges in three dimensions. In particular, we have confirmed to high 
precision that the scaling relation in (|^ holds, and estimated the corresponding amplitude ratio K which we 
believe to be universal. 


A Enumerations 


n 

TAWs 

Bridges 

0 

1 

0 

1 

5 

1 

2 

21 

5 

3 

93 

21 

4 

409 

89 

5 

1853 

369 

6 

8333 

1553 

7 

37965 

6573 

8 

172265 

28197 

9 

787557 

122093 

10 

3593465 

533369 

11 

1647784 

2345429 

12 

7548110 

10366677 

13 

346960613 

46013585 

14 

1593924045 

204927833 

15 

7341070889 

915448621 

16 

33798930541 

4100092693 

17 

155915787353 

18407472565 

18 

719101961769 

82815889677 

19 

3321659652529 

373321398437 

20 

15341586477457 

1685838489629 

21 

70944927549085 

7625255897889 

22 

328054694768261 

34541044018277 

23 

1518490945278377 

156678876463321 

24 

7028570356547189 

711593257794069 

25 

32560476643826933 

3235634079777801 

26 

27 

28 

150838831585499069 

14728414578753489 

67110197685388181 

306074586987649389 


Table 6: Simple cubic lattice TAWs and bridges of length n. 
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Table 7: Simple eubie lattiee irredueible bridges of length n and various heights. 




B Monte Carlo data 


n 

{Xh)n 

{xb)n„ 

511 

0.06923625(26) 

0.1128060(11) 

723 

0.05870409(26) 

0.0956331(11) 

1023 

0.04975890(24) 

0.08104698(73) 

1447 

0.04217388(24) 

0.0686814(10) 

2047 

0.03573665(22) 

0.05818965(67) 

2895 

0.03028231(23) 

0.04930093(93) 

4095 

0.02565517(21) 

0.04176450(63) 

5791 

0.02173570(21) 

0.03538056(85) 

8191 

0.01841284(19) 

0.02996852(60) 

16383 

0.01321296(18) 

0.02150202(57) 

32767 

0.00948020(16) 

0.01542623(53) 

65535 

0.00680166(15) 

0.01106585(50) 

131071 

0.00487920(13) 

0.00793876(47) 

262143 

0.00350038(12) 

0.00569498(44) 

524287 

0.00251118(12) 

0.00408485(40) 

1048575 

0.001801345(94) 

0.00293075(36) 

2097151 

0.001292118(84) 

0.00210227(34) 

4194303 

0.000926791(74) 

0.00150794(30) 

8388607 

0.000665041(66) 

0.00108145(26) 

16777215 

0.000477001(60) 

0.00077565(24) 

33554431 

0.000342261(50) 

0.00055642(23) 


Table 8 : Estimates of {xh)n and {xh)'Hn from Monte Carlo computer experiments. 
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